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In this letter we present the results of our numerical simulations for the finite temperature elec- 
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parameter a s /at ~ 4. The measured thermodynamical quantities (interface tension, jump of 
the order parameter and latent heat) suggest that the phase transition is of very weakly first 
order. 
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1 Introduction 



The standard picture of the electroweak theory tells us that at high temperatures (in the early 
universe) the electroweak symmetry is restored. As the universe expands and cools down, there is 
a phase transition between the high temperature "symmetric" and the low temperature "broken" 
phases. The rate of the baryon violating sphaleron processes is unsuppressed at high temperatures 
and is basically frozen in the low temperature phase. As a consequence, the presently observed 
cosmological baryon asymmetry has been finally determined at the electroweak phase transition [3|. 

In order to clarify the details of this phase transition several techniques have been used. Since 
the bosonic fields have bad infrared behaviour in the perturbative approach 0, 0, || , several nonper- 
turbative approaches have been applied to solve the problem. The most promising one is the use of 
lattice simulations. Two main strategies are used to analyze the problem on the lattice. 

One of them is to use four-dimensional lattices and study the phase transition and its thermody- 
namical properties there. Since the bad infrared behaviour is connected to the bosonic fields only, 
the SU(2)-Higgs model is analyzed || |6|, and the fermionic sector is included by perturbative steps. 
It has been shown that the finite temperature electroweak phase transition is of first order for Higgs 
boson masses around and below 50 GeV. The strength of the phase transition rapidly decreases as 
ran increases. 

The other possibility contains a systematic combination of perturbative and non-perturbative 
methods. One starts with the original theory (e.g. the Standard Model) and integrates out the 
heavy degrees of freedom perturbatively. The obtained theory is a three-dimensional bosonic one 
(e.g. SU(2)-Higgs or SU(2) xU(l)-Higgs) %. 

Analyzing the finite temperature electroweak phase transition by these two different approaches, 
thus in four and three dimensions, provides not only a useful cross-check between them, but the 
comparison of the results also reveals a lot about the applicability of the perturbative reduction 
techniques. Since none of the lattice results contains the fermions, the fermionic sector has to 
be included perturbatively. This would be completely analogous to the perturbative dimensional 
reduction step for the heavy bosonic modes. 

There has been a lot of speculations that the first order nature of the electroweak phase transition 
disappears for large Higgs boson masses (about 80-90 GeV). According to L t = 2, four-dimensional 
results with finite-size analysis there is a first order phase transition for Higgs boson masses of 66 
GeV, but no sign of a first order phase transition could be observed for masses larger than 85 GeV 
(last paper in ||). Recently, in the three-dimensional theory the endpoint of the phase transition 
has been determined ||. Using the same gauge coupling as in the four-dimensional approaches, no 
first order phase transition exists for Higgs boson masses above ~65 GeV. 

It is of essential importance to clarify the relationship between the full four-dimensional results 
and the reduced three-dimensional ones. For this purpose the analysis of the order of the phase 
transition between 65 GeV and 85 GeV Higgs boson masses in the former case is a particularly 
sensitive possibility. The expected shift of the endpoint compared to three dimensions may give 
informations about the overall uncertainties of the reduction procedure. 

In this letter we present the results of our numerical simulations for the finite temperature elec- 
troweak phase transition using the SU(2)-Higgs model on four-dimensional lattices at m# « 80 GeV. 
The temporal extension L t = 2 is used for asymmetric lattice spacings with an asymmetry parameter 
a s /at ~ 4. The measured thermodynamical quantities (interface tension, jump of the order param- 
eter and latent heat) suggest that the phase transition is of very weakly first order. Section 2 gives 
our definitions and some details of the simulations. Section 3 contains the analysis of the thermody- 
namical quantities, while Section 4 presents an estimate of the endpoint for the finite temperature 
electroweak phase transition. In Section 5 we summarize our results and give an outlook. 
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2 Lattice formulation and simulation 



We will study the four-dimensional SU(2)-Higgs lattice model for L t = 2 temporal extensions. Since 
the phase transition is supposed to be quite weak, we need large physical volumes in the simulations. 
The temporal lattice spacing in physical units is set by the critical temperature of the phase transition, 
a t = l/(T c L t ). A large physical volume is ensured by a different choice for the spatial lattice spacing 
a s > a t . The asymmetric lattice spacing scenario can be induced by choosing different coupling 
strengths in the action for time-like and space-like directions. The action reads 

SpM = AE(l-jTrd+A^(l-jTr[/ tp 

sp v L ' tp ^ Z 



+ E^Tr(^)+A 



X 



1 , I 2 
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-k 8 J2 Tr {vt+pU x ^ <p x ) - K t Tr (<p+ +l U xA <p x ) \ , (1) 

where U x ^ denotes the SU(2) gauge link variable, U sp and U tp the path-ordered product of the four 
U x ^ around a space-space or space-time plaquette, respectively; (p x stands for the Higgs field. It is 
useful to introduce the hopping parameter and p 2 = PsPt- The anisotropies 7^ = Pt/Ps and 

7 2 = Kt/n a are functions of the asymmetry £. These functions have been determined with identical 
results on the one-loop level perturbatively || and non-perturbatively JTU] demanding the restoration 



of the rotational symmetry in different physical channels. In this paper we use the specific asymmetry 
parameter £ = 4.052, which gives 7 K = 4 and = 3.919. Details of the simulation techniques can 
be found in ||. 

We choose the parameters as follows. The gauge coupling is set close to its physical value via 
P = 8. As it will be discussed later, its renormalized value receives about 15% correction. The 
scalar self-coupling parameter was chosen to be A = 1.92 • 10~ 4 . For the determination of the critical 
hopping parameter k c and the interface tension (see next section) we have used the two-coupling 
method. The final value for the transition point from this procedure on a 2 x 24 2 x 192 lattice gives 
k c = 0.107791(3). To study final volume behaviour we also simulated on a larger lattice, namely 
2 x 32 2 x 288, where we obtained k c = 0.1077835(25). 

In order to fix the physical parameters in a numerical simulation, one has to define and compute 
suitable renormalized quantities at zero temperature. The renormalized gauge coupling can be de- 
termined from the static potential of an external SU(2) charge pair, measured by Wilson loops. The 
physical Higgs mass mn can be extracted from correlation functions of e.g. R x = |Tr (ip^ip x ) and 
the y9-link operators L^ x = k s T^=i Tr (<f x+jx U X!fl ip x ) + n t Tr {^ +l U xA tp x ) . The W-boson mass m w 
can be obtained similarly from the composite link fields (tp x = p x ■ a x where a x is an SU(2) ma- 
trix and p x > 0) W X]r k = | Tr (r r o& ^U^k a x ^j , r r : Pauli matrices , r, k — 1, 2, 3 . The connected 
correlation functions of these operators have been measured in direction of largest extension on two 
lattices. Note that in principle all numerical estimates for the lightest states in a given channel are 
only upper bounds. Since at our parameters we are deep in the Higgs phase when simulating the 
T = quantities, the lowest mass is well separated from the vacuum and a cross correlation matrix 
analysis does not seem to be very essential. The mass fitting was done by the Michael-McKerrel 



method [11], whose features and application in the SU(2)-Higgs model have been sketched in the 



fourth reference of [|J . The details of the correlation function and Wilson loop analysis for anisotropic 



lattice actions can be found in 10 



At T = and k = k c we obtained the following numerical results on 16 x 8 2 x 64 lattices 
(in units of a t ): the Higgs boson mass is mn = 0.2711(51) and the W-mass is mw = 0.2803(39). 
Combining them with additional data from simulations in the vicinity of k c in order to account for 
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its uncertainty as well, one gets Rhw = rn H /mw = 0.975(50), which corresponds to m H = 78(4) 
GeV pole mass, if mw = 80 GeV sets the physical scale. The renormalized gauge coupling comes 
out to be g\ = 0.539(16). Previous experience shows that the masses and coupling we obtained 
on the above lattice basically coincide with the infinite volume values within our errors. 



3 Thermo dynamical quantities and results 

Near the endpoint the electroweak phase transition falls into the three dimensional Ising universality 
class ||12|| , where the two peak structure persists at any finite volume (even at the endpoint). There- 
fore, we have taken particular attention to study finite volume effects and extrapolate to the infinite 
volume limit. In view of the large lattices one is confronted with for an extremely weak first order 
phase transition also in the case of anisotropic lattices, the interface tension a has been calculated 



by employing the two-coupling method |13] in k. This turned out to be robust and most economic 
in our previous analyses as compared to the histogram or the tunneling correlation length methods 
[[J. Following refs. || 0, the generalization to a situation with anisotropic lattices is straightfor- 
ward. Namely, if an interface pair perpendicular to the ^-direction is enforced by dividing the lattice 
volume in symmetric and Higgs phases with 

k = k 2 ) = (ki < k c : z < L t /2 , k 2 > k c : z > L t /2) (2) 

as for the ^-determination, the related additional free energy AF yields for An = \Ki — K c \ <C 1 the 
expression 

a 2 s a t a = lim <^ (An) lim L z ■ AL v (ki, k 2 ) > , (3) 



Are— >0 I L z — >oo 



where = L^(ki, k 2 ), i = 1,2, denote the expectation value of the <£>-link operators in the 
respective phases and AL v (ki, k 2 ) = L^(«i, n 2 ) — L^(ki, k 2 ) their difference. Now, since the free 
energy shift can be shown to behave as AF w 0{Ak), the (N + 2)-parametric Laurent ansatze 

N 

L^(K U K 2 ) = --^ + Y,li J \^-K c y + 0((AK) N + 1 ) ,2 = 1,2 (4) 

give the finite volume estimator for the interface tension a 2 s a t a = L z (ci + c 2 ) . 

In the left part of figure [1] we illustrate a characteristic two-phase distribution of ^-slices for the 
expectation value L v (z). An interface pair has developed, and the plateaus marking the pure-phase 
expectation values L^(ki,k 2 ) are flat and broad enough to ensure that the coexisting phases are 
stable against any turn-over into one single phase. All simulation data, which were collected in 
k = k 2 ) centred symmetrically around k c for the determination of a, are listed in table |l|. The 
integrated autocorrelation times of the <£>-links are 18 to 160 sweeps for our largest and smallest 
Ac-interval, respectively. Compared to previous analyses on isotropic lattices for stronger first order 
transitions, the ^-intervals have been chosen closer to each other and to k c , while it was always 
verified that the distinct phases were well separated. At a very weak transition the diverging of L$ , 
when K\ (k 2 ) approach k c from below (above), sets in later and will be less pronounced than for 
smaller Higgs masses. Therefore, the curvatures in the shapes of the two 1$ as functions of k and 
particularly their residua should be modelled more accurately in the present case. The least-squares 
fits to eq. (|j) can be done for both if- links individually or for their difference AL 9 . We observed || 
that the second alternative should be preferred. Firstly, correlations between the competing phases 
give fluctuations in the locations of the interfaces and secondly, a roughening of the interfaces is 
expected as the phase transition weakens. The two effects are supposed to compensate to some 
extent in AL V . 
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40000 
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0.107771 


0.107811 


60000 


13.1186(89) 


20.7117(90) 


7.593(11) 


0.107776 


0.107806 


60000 


13.3410(90) 


20.215(13) 


6.874(12) 


0.107781 


0.107801 


70000 


13.711(15) 


19.599(16) 


5.888(20) 


0.107786 


0.107796 


80000 


14.316(39) 


18.773(25) 


4.457(37) 



Table 1: Results for Lffi and AL V for the calculation of a on a lattice of size 2 x 24 2 x 192, 
together with their statistical errors from a binning procedure. In order to obtain acceptable x 2 - 
values, the entries in italics were only omitted for the four-parameter fits (see text). 



In this spirit we carefully examined nearly all types of fits from three- to six-parametric in the 
Laurent ansatze of eq. @j for the two kinds of observables and AL V and for every allowed fit 
intervals, which could be built up. In part the results on a might depend on the specific fit interval 
and on the number of fit parameters in question, but if the inverse-linear part of the general ansatz 
is supplemented with a fourth or even a fifth fit parameter, the fits are acceptable throughout with 
tolerable x 2 /dof and can readily be extended over almost the whole range of K-pairs. In order to 
control the relevance of a further degree of freedom, we also performed some fits with an additional 
parameter multiplying the next higher power in — k c . This gave either unreasonably too large 
X 2 /dof- values, or the error analysis of these fits led to very high statistical errors so that the new fit 
parameters could not be resolved reliably. As exemplarily displayed in the right diagrams of figure [I] 
and in figure || we found for the best fits (a/Tj}) L ^ = 0.00060(8 + 22) from four-parameter x 2 -fits 

of L«, i = 1,2, and (<r/T c 3 ) AL ^ = 0.00060(7) from the similar fit of AL V . The quoted error on 
the former incorporates the statistical error from a bootstrap analysis |T4] , |6| and an error caused by 
the uncertainty in k c (second entry in parentheses), whereas the solely statistical error on the latter 
stems from repeated fits with 1000 normally distributed random data around the measured mean 
values as input. 

To study possible finite volume effects and the infinite volume limit, we also simulated on a larger 
lattice of size 2 x 32 2 x 288 at otherwise unchanged parameters and couplings. The data are displayed 
in table 2. After an analogous evaluation we obtained for (a/T^) L ^ = 0.00062(10 + 25) from four- 
parameter x 2 ~fits of i = 1,2, and (d'/T^) AL = 0.00063(8) from the similar fit of AL V . These 
are fully compatible within errors to the preceding results, which indicates that the surface tension 
does not decrease for increasing volume when the volume is encreased by a factor of 2.7. Thus data 
are compatible with a first order phase transition. 

To present our final (infinite volume) value we average all these results to a combined estimate and 
take their absolute spread including the errors as a measure for the sum of statistical and systematic 
uncertainties. Assuming all reasonable fits from both spatial volumes to contribute, we finally arrive 
at 

all 

= 0.0006(4) . (5) 

1-K 

From quadratic fits of the discontinuities of the order parameters showing up in the thermal cycles 
of figure for p 2 and L v (and similarly for the plaquette variables) we also extracted the jump in 
the Higgs field vacuum expectation value and the latent heat. The jump of the Higgs field vacuum 




expectation value is given by Av/T c = L t £ V2kA(p 2 ) in case of asymmetric lattice spacings. The 
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Figure 1: Left: two-phase profile of the z-slice expectation value L v {z) of L v for the smallest 
Ac-interval used in the simulations for the a-determination in table |l| Right: four-parameter least- 
squares fits of L®, i = 1,2, separately in each phase. The x 2 ~ va hies are xi — 2.59 and x\ — 0.96, 
respectively. 



latent heat is determined by 



Ae 

y4 



L\ OAs 



(6) 



where A means the jump, s is the average action density and r = — ln(a s mw) and the derivative 
is taken with constant £. To determine the partial derivatives of /3 S , (3 t and A we have used the 
perturbative renormalization group equations defining the lines of constant physics, while for k we 
have followed a procedure similar to the one used in [Hj for the isotropic case. We end up with 



Ae 

c 



1 dn A , 43 
-— A < i v > + 
k or 



16/T 2 



<? 2 (&A < Ps > +/3 t A <p t >) 



. „ (2\3k 
+A < Q > — — + 



K <9t 4^7T 2 



96A C 2 + - 9A C # 2 



(7) 



where < > is the average link, < p s > and < p t > are the average plaquette variables, < Q > is 
the average of the quartic part in (p of the lattice action (|l|) and A c = A£/(4k 2 ). The results for our 
simulation parameters are collected in table |^. 



4 Estimate for the endpoint for the electroweak phase tran- 
sition 

It is a well known perturbative feature that the finite temperature electroweak phase transition is 
of first order due to the vanishing magnetic mass (%). A non- vanishing parametrization of this 
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Figure 2: Four-parameter least-squares fit of the <y9-link difference AL V as a function of k-i — K\ 
with x 2 — 4.55. The extrapolation to k% — K\ — gives the interface tension. 



mass predicts an endpoint mn, c f° r the first order phase transition. Note that our treatment is 
purely phenomenological and is performed in one particular gauge, the Landau gauge. The value 
of the magnetic mass can be obtained e.g. by solving the truncated Dyson-Schwinger equations 
(gap-equations) [|3], [HJ. Different treatments of these equations give different numerical values for 



the magnetic mass. Nevertheless keeping it as an unknown parameter and analysing the behaviour 
of the different static thermodynamical quantities (e.g. jump of the order parameter, latent heat 
and interface tension), a scaling behaviour can be observed in the vicinity of vnn c ||. As we have 
explicitly checked, this behaviour is essentially independent of the model or loop-order used |4j]. 

We have fitted all our available four-dimensional data at m H « 19, 35, 49, 80 GeV (see |J and the 
results of the present paper) on the static thermodynamical quantities to the predictions of the one- 
and two-loop finite temperature effective potential. The perturbative results were considered as a 
function of one unknown parameter, namely the endpoint of the phase transition m# jC parametrized 
by a non-vanishing mass in the magnetic sector. The small masses play a less relevant role in the 
fiting procedure, while the influence of 49 and 80 GeV data on the value of the endpoint is more 
important. Combining all results (i.e. all masses and thermodynamical quantities), we have seen that 
at the optimal magnetic mass x 2 /dof ~ 1, whereas for vanishing mj- the x 2 -value is considerably 
larger. This fact can be interpreted as a sign, which indicates the presence of the phase transition 
endpoint. Our combined value for the endpoint is m# >c = 106l|f GeV. The estimated value for 
m# )C is slightly larger than the values of the three-dimensional analyses ||; however, they are in an 
approximately 1.5— a agreement. We emphasize that at small L t 's the phase transition usually tends 
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Table 2: Results for Lffi and AL V for the calculation of a on a lattice of size 2 x 32 2 x 288, 
together with their statistical errors from a binning procedure. 





T c /m H 




&v/T c 


A e /T c 4 
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1.86(2) 


0.0006(4) 


0.37(16) 


0.0033(27) 



Table 3: Summary of our lattice results at L t = 2 temporal extension. 



to be stronger than for larger L t values . 

5 Discussion and outlook 

Both quantities cr/T c 3 and Ae/T c 4 for itlh ~ 80 GeV are substantially smaller than the perturbative 
predictions with zero magnetic mass (e.g. cr/T 3 « 0.002 [0). Comparing with our earlier investi- 
gations at lower itlh f| §], these results confirm the expectation that the interface tension and the 
latent heat are steeply decreasing functions of the Higgs boson mass, and they are even consistent 
at mu ~ 80 GeV with a no first order phase transition scenario approximately at the 1—a level. 

We have fitted all our four-dimensional data to the perturbative predictions assuming, however, 
an endpoint for the phase transition parametrized by a phenomenological non-vanishing magnetic 
mass. Our combined value for the endpoint is m^ c = 106^24 GeV. Although being in a nearly 1.5-er 
agreement, the fact that our findings deviate from those of the three-dimensional investigations 0, 
which claim the endpoint of the first order phase transition line to be at m# )C ~ 67 GeV, should be 
clarified in future. 

However, we also have to concede that a temporal lattice extension of L t = 2 may be still too 
far from the continuum physics. In fact, experience with lattice perturbation theory shows |T(J that 



using anisotropic lattices makes the approach to the continuum limit even slower than for isotropic 
lattices. So at least some knowledge about the behaviour of the four-dimensional SU(2)-Higgs model 
at run ~ 80 and L t = 3 is of principal interest before drawing a final conclusion. Unfortunately 
the necessary CPU-time requirements to reach an adequate precision on such huge lattices seems 
not yet to be realistic. Therefore, it would be more desirable to determine the endpoint of the 
electroweak phase transition on the basis of a study of the Lee- Yang zeros, as applied in refs. || to 
the dimensionally reduced model, within the theory in four dimensions as well. 
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